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1.  Introduction 


The  pitch-damping  moment  coefficients  Cmq  (due  to  body  transverse  angular  rate)  and  C 

(due  to  angular  rate  associated  with  angle  of  attack)  play  an  important  role  in  the  performance 
and  dynamic  stability  of  flight  bodies.  The  pitch-damping  moment  coefficient  sum  Cmq  +  Cmd, 

is  of  most  practical  importance,  although  the  individual  damping  coefficients  are  often  required 
in  aerodynamic  analyses.  Throughout  the  past  several  decades,  a  variety  of  techniques  and 
theories  have  been  developed  for  predicting  the  pitch-damping  coefficients  (1-10).  These 
techniques  vary  in  their  ease  of  use  as  well  as  their  ability  to  accurately  predict  the  pitch¬ 
damping  coefficients. 

During  the  course  of  its  development,  slender  body  theory  was  generalized  to  predict  a  large 
variety  of  aerodynamic  coefficients  including  the  pitch-damping  coefficients  (1-2).  In  general, 
direct  application  of  these  methods  provides  only  qualitative  results  for  the  aerodynamic 
coefficients.  However,  elements  of  slender  body  theory  have  been  incorporated  into  current 
engineering  methods.  These  methods  (3-5)  have  evolved  considerably,  although  their 
implementation  is  fairly  complex.  Apart  from  implementation  issues,  modern  engineering 
methods,  once  embodied  into  a  computer  code,  are  relatively  easy  to  use  and  provide  fast  and 
reasonably  accurate  aerodynamic  predictions  for  a  large  variety  of  flight  geometries. 


From  slender  body  theory,  some  important  relationships  between  the  various  aerodynamic 
coefficients  can  be  derived,  although  these  relationships  only  hold  rigorously  within  the  context 
of  theories  from  which  they  were  obtained.  Bryson  (2)  derived  the  relatively  well-known  slender 
body  result  that  relates  the  pitch-damping  moment  coefficient  sum  to  the  normal  force 
coefficient,  shown  in  equation  1 : 
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(1) 


Sacks  (1 ),  using  the  Blasius  method  for  calculating  the  forces  and  moments  on  slender  bodies 
from  the  cross-flow  potential,  found  that  many  of  the  aerodynamic  coefficients  were  related  to 
each  other.  Sacks  obtained  expressions  that  directly  related  the  individual  pitch-damping 
coefficients  to  the  nonnal  force  coefficient,  such  as  the  pitch-damping  force  coefficient  shown  in 
equation  2.  These  expressions  have  a  form  similar  to  Bryson’s  result  shown  in  equation  1 : 


C 


Nn 


-x  A 

ACg 

D 


C 


(2) 


In  practice,  these  relations  that  directly  relate  the  damping  coefficients  to  the  normal  force 
coefficient  do  not  perform  particularly  well,  even  when  the  slender  body  evaluation  of  the 
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normal  force  coefficient  is  replaced  with  a  more  accurate  evaluation  of  the  normal  force 
coefficient  from  sources  such  as  experimental  measurement  or  computational  fluid  dynamics 
(CFD)  (6).  However,  these  relations  can  be  combined  with  empirical  corrections  to  yield  more 
reliable  results  (6,  7). 

Sacks  also  found  expressions  that  related  the  individual  damping  coefficients  to  each  other, 
including  the  following  relationship  between  the  pitch-damping  force  coefficients: 

CNq=CNa-Cma.  (3) 


Sacks’  explicitly  derived  relation  shown  in  equation  3  can  be  easily  generalized  using  his  theory 
to  the  individual  pitch-damping  force  and  moment  coefficients,  and  the  pitch-damping  force  and 
moment  sums  as  shown  in  equations  4-9.  For  the  purposes  of  this  report,  these  relationships  will 
be  referred  to  as  Sacks’  relations: 
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Sacks’  contribution  is  the  recognition  that  these  relationships  exist,  although  the  same  relations 
are  implicitly  contained  in  other  theories  such  as  that  presented  by  Bryson  (2). 


For  slender  bodies,  these  relations  are  independent  of  configuration  and  are  applicable  to  both 
winged  and  wingless  bodies.  The  importance  of  these  relationships  is  that  if  just  one  of  the  three 


pitch-damping  force  (or  moment)  coefficients  can  be  determined  ( CN  ,  CN  or 
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the  other  two  damping  coefficients  can  be  obtained  using  simple  closed-form  expressions.  This, 
of  course,  assumes  that  the  pitching  moment  coefficient  Cma  (first  moment  of  the  normal  force) 

and  the  second  moment  of  normal  force  ^  2ma  can  be  obtained  as  well.  Both  of  these 
coefficients  can  be  obtained  if  the  normal  force  distribution  is  known  as  shown  in  equations  10 
and  1 1 : 
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Predictive  methods  for  the  static  normal  force  and  static  pitching  moment  are  well  established, 
even  for  fast-design  methods. 

Because  Sacks’  relations  were  derived  using  simple  approximate  theories,  it  remains  to  be  shown 
whether  the  validity  of  the  relations  shown  in  equations  4-9  exists  only  within  the  context  of  the 
theories  from  which  they  were  derived  or  whether  they  are  universally  valid  for  slender  bodies, 
or  perhaps,  more  importantly,  whether  they  are  of  general  engineering  significance.  It  is 
important  to  note  that  the  relations  shown  in  equations  4-9  differ  somewhat  from  the  results  in 
equations  1  and  2  because  they  relate  the  pitch-damping  coefficients  to  each  other  rather  than 
purely  to  the  static  normal  force  coefficient.  This  suggests  that  Sacks’  relations  may  more 
properly  represent  the  physics  compared  to  their  counterparts  which  directly  relate  the  damping 
coefficients  to  the  nonnal  force  coefficient.  Because  of  this,  there  is  reason  to  investigate  the 
validity  and  accuracy  of  Sacks’  relations. 

Recently,  a  computational  approach  for  predicting  all  three  of  the  pitch-damping  coefficients  has 
been  developed  (5,  9).  The  approach  solves  the  three-dimensional  (3-D)  thin-layer  Navier- 
Stokes  equations  for  three  different  imposed  motions  that  allow  the  three  pitch-damping 
coefficients  to  be  predicted  independently.  The  pitch-damping  force  and  moment  coefficient 
sums  are  detennined  from  the  computation  of  a  body  undergoing  an  imposed  coning  motion. 

The  individual  pitch-damping  coefficients  are  obtained  from  computations  of  a  body  undergoing 
two  specific  types  of  imposed  helical  motions.  Each  of  these  motions  is  described  in  the 
following  sections.  One  of  the  key  components  of  this  method  is  that  steady  flow  techniques  can 
be  employed  to  predict  aerodynamic  derivatives  normally  associated  with  time-dependent 
motions. 


1.1  Helical  Motions 

Forces  and  moments  related  to  the  two  individual  rates  q  and  a  can  be  excited  independently 
using  two  types  of  motion  where  the  center  of  gravity  (CG)  of  the  flight  vehicle  traverses  a 
helical  flight  path.  The  first  motion  requires  the  vehicle’s  longitudinal  axis  to  be  oriented  in  the 
same  direction  as  the  center  of  rotation  of  the  helix  but  displaced  by  a  constant  distance. 

Figure  1  shows  a  3-D  view  of  the  motion. 

This  particular  motion  produces  no  rotation  of  the  body-fixed  nonrolling  coordinate  frame 
relative  to  an  earth-fixed  coordinate  frame,  and  hence,  the  transverse  angular  velocity  of  the  body 
is  zero.  The  angle  of  attack  and  its  angular  rate  vary  continuously,  producing  moment 
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components  associated  with  the  coefficients  Cm^  and  Cm^  ,  respectively.  This  motion  is 

referred  to  as  “c/  =  0  helical  motion ”  because  the  angular  rates  associated  with  the  damping 
coefficient  Cm  are  zero. 

For  the  second  motion,  the  longitudinal  axis  of  the  flight  vehicle  remains  tangent  to  the  helical 
flight  path  at  each  point  along  the  trajectory.  Figure  2  shows  a  3-D  view  of  this  motion.  The 
angle  of  attack  of  the  incident  airstream  is  zero  because  both  the  longitudinal  axis  of  the  body 
and  the  free-stream  velocity  vector  are  tangent  to  the  flight  path.  The  resulting  yawing  rate  is 
also  zero  because  the  angle  of  attack  is  constant.  The  angular  orientation  of  the  flight  body 
changes  continuously  with  respect  to  an  earth-fixed  reference  frame,  producing  a  nonzero 
transverse  angular  rate. 

As  a  result,  moment  components  associated  with  the  damping  moment  Cmq  are  produced.  This 
motion  is  referred  to  as  “  a  =  0  helical  motion ”  because  the  angular  rates  associated  with  the 
damping  coefficient  Cm^  are  zero. 

For  each  of  the  helical  motions,  the  transverse  aerodynamic  moment  in  the  nonrolling  frame  will 
be  periodic  in  time,  which  also  indicates  that  the  flowfield  will  be  periodic  in  time  when  viewed 
from  the  nonrolling  coordinate  frame.  The  time  dependency  is  removed  by  transfonning  to  an 
orthogonal  right-handed  coordinate  system  that  has  its  x-axis  aligned  with  the  longitudinal  axis 
of  the  body  and  its  z-axis  along  a  line  between  the  body’s  CG  and  the  axis  of  rotation  of  the 
helix. 
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Figure  2.  Helical  motion  with  zero  <2  and  nonzero  q  . 


For  each  of  the  helical  motions  previously  described,  the  spin  rate  of  the  body  has  not  been 
defined.  To  eliminate  any  contributions  to  the  aerodynamic  forces  and  moments  from  the 
Magnus  forces  and  moments,  the  spin  rate  is  fixed  to  zero  (see  references  [5]  and  [9]  for  details). 
The  resulting  in-plane  moment  (Cm )  and  side  moment  (Cn )  coefficients  in  the  transfonned 
coordinate  system  for  both  type  of  helical  motions  are  shown  in  equations  12  and  13: 

Zero  -  spin  q  =  0  helical  motion 


C  +  i  C  =  -C 
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Zero  -  spina  =  0  helical  motion 


C  +i  C  =C 


QD 
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R0Q 


V 


(13) 


Here,  Q  is  the  angular  velocity  of  the  body  about  the  helix  axis,  R0  is  the  perpendicular 

distance  between  the  helix  axis  and  the  body  CG,  and  V  is  the  total  linear  velocity  of  the  CG. 
Similar  expressions  for  the  individual  damping  force  coefficients  can  be  developed  using  the 
same  approach  as  applied  for  the  moment  coefficients. 


1.2  Coning  Motion 

To  predict  the  pitch-damping  coefficient  sum,  coning  motion  is  employed.  In  steady  coning 
motion,  the  longitudinal  axis  of  the  flight  body  performs  a  rotation  at  a  constant  angular  velocity 
about  a  line  parallel  to  the  free-stream  velocity  vector  and  coincident  with  the  body’s  CG,  while 
oriented  at  a  constant  angle  with  respect  to  the  free-stream  velocity  vector.  This  is  shown 
schematically  in  figure  3.  In  the  context  of  this  report,  coning  motion  also  requires  the  CG  to 
traverse  a  rectilinear  path  at  constant  velocity  such  that  the  free-stream  velocity  vector  has  a 
fixed  orientation  with  the  inertial  frame. 


Relative  to  a  nonrolling  coordinate  frame,  both  the  transverse  angular  rate  of  the  body  and  the 
angular  rate  associated  with  the  angle  of  attack  vary  in  a  periodic  manner  thereby  exciting  the 
aerodynamic  forces  and  moments  associated  with  both  of  the  individual  pitch-damping 
coefficients.  Here,  a  specific  form  of  coning  motion,  described  as  zero-spin  coning  motion,  is 
employed,  in  zero-spin  coning  motion,  the  total  angular  velocity  of  the  body  along  the 
longitudinal  axis  (the  spin  rate)  is  zero.  By  imposing  zero  spin  rate  on  the  body,  the 
contributions  from  the  Magnus  forces  and  moments  are  eliminated. 

The  time  dependency  is  removed  by  transforming  the  body-fixed  nonrolling  coordinate  frame  to 
an  orthogonal  right-handed  coordinate  system  that  has  its  x-axis  aligned  with  the  longitudinal 
axis  of  the  body  and  its  z-axis  in  the  pitch-plane  of  the  body.  Within  this  transfonned  coordinate 
frame,  in-plane  moment  (Cm )  and  side  moment  (Cn )  coefficients  have  the  following  form: 
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Zero-spin  coning  motion 
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Here,  the  side-moment  is  proportional  to  the  pitch-damping  moment  coefficient  sum  and  varies 
linearly  with  the  coning  rate  Q  and  sine  of  the  total  angle  of  attack  8.  For  small  angles  of 
attack,  the  cosine  of  the  total  angle  of  attack,  y,  can  be  assumed  to  be  one. 


The  pitch-damping  coefficient  sum  can  also  be  determined  by  simply  adding  the  individual 
damping  coefficients.  In  practice,  there  is  very  little  difference  in  the  two  results.  However,  for 
the  current  study,  directly  predicting  the  pitch-damping  coefficient  sum  using  coning  motion 
provides  an  alternative  prediction  of  the  pitch-damping  sum  and  additional  confirmation  of  the 
predictions  of  the  individual  coefficients. 


2.  Computational  Technique 


In  the  previous  section,  several  types  of  steady  motion  were  presented  that  produce  aerodynamic 
forces  and  moments  from  which  the  various  pitch-damping  coefficients  can  be  obtained.  One 
unique  feature  of  these  motions  is  that  they  are  steady.  The  advantage  of  a  steady  motion  over 
an  unsteady  motion  is  that  a  potentially  time-independent  flowfield  can  be  produced  by  a  steady 
motion,  permitting  analysis  using  steady  flow  CFD  techniques.  Such  techniques  can  be 
computationally  less  expensive  than  time-dependent  CFD  approaches.  To  fully  exploit  the 
steady  character  of  the  flow,  special  body-fixed  coordinate  systems  have  been  employed  to 
capture  the  steady  flowfield.  One  feature  of  these  coordinate  frames  is  that  they  are  rotating  at  a 
constant  rate  with  respect  to  an  inertial  frame.  Because  of  this,  the  governing  equations  of  fluid 
motion  must  be  modified  to  take  into  account  the  centrifugal  and  Coriolis  force  terms  associated 
with  the  noninertial  rotating  frame. 


The  steady  thin-layer  Navier-Stokes  equations  are  shown  in  equation  15: 


GF  GF  GG  -  -  1 

—  +  —  + +  Hn  +  H  =  — 

GE,  drj  dC,  Re 


*+s. 

ec  0 


(15) 


The  inviscid  flux  vectors,  E ,  F ,  and  G  ,  the  viscous  term  S  ,  the  inviscid  and  viscous  source 
terms  due  to  the  cylindrical  coordinate  formulation,  Hc  and  Sc ,  and  the  source  term  H, 
containing  the  Coriolis  and  centrifugal  force  terms  that  result  from  the  rotating  coordinate  frame, 
are  functions  of  the  dependent  variables  represented  by  the  vector  qT  =  (p,  pu,  pw,  pm,o) ,  where 
p  and  e  are  the  density  and  the  total  energy  per  unit  volume,  and  u ,  v ,  and  w  are  the  velocity 
components  in  axial,  circumferential,  and  normal  directions.  The  inviscid  flux  vectors  and  the 
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source  term  are  shown  in  equation  16.  Details  of  the  thin-layer  viscous  term  are  available  in  the 
literature  (11): 


pU 

pV 

PuU  +  £xp 

1 

F  =  - 
J 

puV  +  ?7xp 

pv  U 

pvV  +  rj^/r 

pwU 

pwV  +  rjrp 

(e  +  p)U  _ 

(e  +  p)V 

pW 

0 

PuW  +  cxp 

H  =  t 
J 

pU 

PvW  +  ^p/r 

pU 

pwW  +  <^rp 

pfr 

(e  +  p)W 

_/?ufx  +/?vf^  +/?wfr_ 

(16) 


The  pressure,  p,  can  be  related  to  the  dependent  variables  by  applying  the  ideal  gas  law: 


p  =  (r-i) 


;_£(u2+v2+w2' 


(17) 


The  turbulent  viscosity,  pt,  which  appears  in  the  viscous  matrices,  was  computed  using  the 
Baldwin-Lomax  turbulence  model  (12). 

The  Coriolis  and  centrifugal  acceleration  terms  due  to  the  rotating  coordinate  system,  which  are 
contained  in  the  source  tenn,  H,  are  shown  in  equation  18: 

f  =  2Q  x  u  +  Q  x  (Qx  r)  .  (18) 


The  Coriolis  acceleration  is  a  function  of  the  angular  velocity  of  the  coordinate  frame  with 
respect  to  the  inertial  frame,  Q,  and  the  fluid  velocity  vector,  u,  which  can  be  represented  by  the 
velocity  components,  u ,  v ,  and  w .  The  centripetal  acceleration  is  a  function  of  the  angular 
velocity  of  the  rotating  frame,  Q,  and  the  displacement  vector,  R,  between  the  axis  of  rotation 
and  the  local  position  in  the  flowfield.  The  acceleration  vector,  f,  can  be  written  in  terms  of  its 
components  along  the  x ,  (/>,  and  r  axes,  fx ,  f(4 ,  and  fr . 

The  steady  thin-layer  equations  are  solved  using  the  parabolized  Navier-Stokes  (PNS)  technique 
of  Schiff  and  Steger  (11).  This  “space-marching”  approach  integrates  the  governing  equations 
from  the  nose  of  the  flight  body  to  the  tail.  Following  the  approach  of  Schiff  and  Steger,  the 
governing  equations,  which  have  been  modified  here  to  include  the  Coriolis  and  centrifugal  force 
terms,  are  solved  using  a  conservative,  approximately  factored,  implicit  finite-difference 
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numerical  algorithm  as  fonnulated  by  Beam  and  Warming  (13).  Details  of  the  implementation 
of  the  source  tenn  that  contains  the  Coriolis  and  centrifugal  force  terms  are  given  in  references 
(8)  and  (9). 

The  technique  has  been  validated  with  available  experimental  data  where  possible,  and  excellent 
agreement  is  found  (8,  14).  Grid  resolution  studies  were  also  performed  in  the  original  studies  to 
ensure  grid-independent  solutions  (8,  9). 


3.  Results 


In  the  current  context,  this  computational  procedure  allows  the  general  applicability  of  Sacks’ 
pitch-damping  relations  to  be  examined.  Arguably,  up  to  now,  it  has  not  been  possible  to  assess 
the  validity  or  accuracy  of  these  relationships  due  to  the  uncertainty  associated  with 
experimentally  derived  pitch-damping  data  and  the  lack  of  a  higher  order  theory.  In  the  current 
research  effort,  the  technique  has  been  applied  to  two  axisymmetric  body  geometries:  an  ogive- 
cylinder  configuration  and  a  cone-cylinder  body  with  a  flared  afterbody.  The  accuracy  of  Sacks’ 
relations  is  examined  for  each  of  these  body  geometries  in  supersonic  flight. 

3.1  Ogive-Cylinder  Results 

The  computational  approach  was  applied  to  a  secant-ogive  cylinder  body  (designated  as  the 
Army  Navy  Spinner  Rocket  [ANSR])  shown  in  figure  4.  Results  for  sea-level  atmospheric 
conditions  ( RcA)  =  4.53  x  105  A/x  )  are  shown  as  an  example,  although  other  Reynolds  numbers 
and  nose  geometries  were  considered  during  the  course  of  the  study.  Baseline  results  are  for  a 
flight  velocity  of  Mach  2.5  with  additional  results  at  Mach  1.8,  3.5,  and  4.5.  The  results 
presented  here  are  representative  of  the  other  flight  conditions  and  nose  geometries  examined. 

For  the  force  coefficients,  the  deviation  or  “error”  in  the  application  of  Sacks’  relations  is  defined 
in  equation  19.  The  deviation  Aj  represents  the  difference  between  the  right-hand  and  left-hand 

sides  of  equation  4: 

A  =  CNa-Cma-CNq  =A,  (19) 


Note  that  by  simple  algebraic  manipulations,  the  following  relations  that  are  algebraically 
equivalent  to  equation  19  are  found: 


and 


A  = 


CNn  +C 


tc 


-Cm  -2CNn  =  A 


k2  ’ 


(20) 


A  =  2Cn  - C 


+  C 


N/- 


(21) 
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Essentially,  these  deviations  are  a  measure  of  the  accuracy  of  the  Sacks’  relations  shown  in 
equations  4,  6,  and  8.  It  could  be  argued  that  the  differences  A,,  A2,  and  A3  are  numerical 
errors  that  are  not  representative  of  physical  phenomena.  However,  it  should  be  noted  that  the 
deviation  A  can  be  computed  independently  in  three  different  ways.  The  deviation  A|  is 
computed  using  the  individual  damping  coefficients  CNq  and  C  obtained  independently 
from  two  different  types  of  helical  motion,  while  A2  and  A3  are  computed  from  the  pitch¬ 
damping  coefficient  sum  CN  +  CN  .  and  one  of  the  two  individual  damping  coefficients  CN 
and  CN  ,  respectively.  The  importance  of  having  three  independent  methods  for  determining  A 
is  to  demonstrate  that  the  deviation  is  not  due  to  an  “error”  in  the  prediction  of  any  one  of  the 
three  damping  coefficients.  Similar  expressions  for  the  deviation  of  the  moment  coefficients  A 
can  also  be  obtained  and  are  shown  in  equations  22-24: 

A  =  Cm.  -C^-C^Aj,  (22) 


and 


A  = 


+  CT 


_C2m  -2Cm_  =  A 


a 


2  ’ 


(23) 


A=2Cm.  -C 


2mr 


+  cr 


(24) 
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Figure  5  shows  the  computed  deviation  for  the  force  coefficient  as  a  function  of  body  length  at 
Mach  2.5.  Here,  the  body  has  been  lengthened  to  20  calibers  to  more  clearly  illustrate  the 
variation  of  the  deviation  with  body  length.  There  is  very  good  correlation  between  the  three 
different  methods  of  computing  A ,  suggesting  that  the  deviation  is  representative  of  a  physical 
effect  rather  than  simply  numerical  error.  The  deviation  is  very  small  over  the  nose  and  begins 
to  grow  in  a  linear  fashion  a  couple  of  body  diameters  aft  of  the  nose.  The  deviation  (for  the 
force)  can  be  shown  to  be  independent  of  CG  position  using  the  CG  translation  relations  for  the 
individual  force  coefficients.  This  was  also  confirmed  by  varying  the  CG  position  in  the 
computations  as  well.  The  computed  deviation  was  also  found  to  be  somewhat  dependent  on 
nose  length  and  Mach  number,  although  the  results  shown  are  representative  of  the  trends 
observed  for  the  other  flight  conditions. 


0  4  8  12  16  20 

X/D 


Figure  5.  Comparison  of  the  three  computed  deviations  At ,  A2 ,  and  A3  vs.  body 
length,  ANSR,  and  Mach  2.5. 

Figure  6  shows  the  computed  deviations  for  the  force,  A,,  A,,  and  A3,  compared  with  the 
individual  pitch-damping  coefficients  and  the  pitch-damping  sum.  Relative  to  CN^  ,  and  the 
pitch-damping  force  sum,  the  deviation  is  quite  small.  The  distribution  of  the  deviation  is  also 
small  compared  with  the  pitch-damping  force  coefficient  CNq  ,  except  near  the  end  of  the  body 
where  CNq  itself  is  nearly  zero.  Equation  4  shows  that  according  to  Sacks’  relation,  the 
difference  between  CN^  and  CNq  is  the  pitching  moment  coefficient  Cm  .  This  is  also 
graphically  shown  in  figure  6.  Near  the  nose  of  the  body,  the  difference  between  the  two 
damping  coefficients,  CN^  and  CNq  ,  increases  at  nearly  the  same  rate  as  the  pitching  moment 
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Figure  6.  Comparison  of  the  three  computed  deviations  A, ,  A2 ,  and  A3  relative  to  the 
damping  force  coefficients,  ANSR,  Mach  2.5,  x  /L  =  0.5598  ,  and  length-to- 
diameter  (L/D)  ratio  =  9. 

coefficient.  On  the  aft  of  the  body,  the  pitching  moment  coefficient  is  nearly  constant,  and  the 
difference  between  CN^  and  CNq  becomes  relatively  constant. 

Figure  7  shows  the  computed  deviation  for  the  moment  relative  to  the  three  damping  coefficients 
of  interest.  Again,  there  is  very  good  correlation  between  the  three  methods  for  computing  the 
deviation.  The  computed  deviation  is  small  compared  to  individual  damping  moment 
coefficients  and  to  the  pitch-damping  moment  sum.  Sacks’  relation  for  the  moment  (equation  5) 
shows  that  the  difference  between  Cm^  and  Cm^  is  the  pitching  second-moment  coefficient 
C2m  .  This  also  can  be  seen  in  figure  7.  In  a  sense,  figures  6  and  7  provide  some  measure  of 
the  expected  error  in  applying  Sacks’  relations  to  obtain  the  various  damping  moment 
coefficients  given  that  normal  force  distribution  (or  Cma  and  C2ma  )  and  one  of  the  pitch¬ 
damping  coefficients  is  known. 

The  most  likely  application  of  Sacks’  relations  in  practical  situations  is  to  compute  the  individual 
damping  coefficients  from  the  pitch-damping  coefficient  sum  because  the  pitch-damping 
moment  sum  is  much  easier  to  measure.  Through  simple  algebraic  manipulations  of  equations 
7  and  9,  the  following  form  of  the  Sacks’  relations  can  be  obtained: 
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(25) 


and 


+  Cr 


+  C 
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+  Cr 


(26) 


Figure  7.  Comparison  of  the  three  computed  deviations  A, ,  A2 ,  and  A3  relative  to  the 

damping  moment  coefficients,  ANSR,  Mach  2.5,  xcg/L  =  0.5598 ,  and  L/D  =  9. 

Figure  8  shows  a  comparison  of  the  pitch-damping  moment  coefficients  Cm^  and  Cm^ 

obtained  by  applying  the  Sacks’  relationships  using  the  pitch-damping  moment  coefficient  sum 
and  the  second-moment  coefficient  of  the  normal  force  C2m  .  Here,  CFD  has 


Cmq  +  C  m  (y 


been  used  to  compute  both 


Cmq  +  Cma 


and  C2m  .  The  results  from  the  application  of  these 


relations  are  compared  with  the  CFD  predictions  of  the  individual  damping  coefficients.  The 
predictions  were  obtained  at  a  flight  velocity  of  Mach  2.5.  The  coefficients  Cm  and  Cm^  are 

overpredicted  and  underpredicted  by  ~5%  and  12%,  respectively,  although  the  absolute  error  is 
similar  for  both  coefficients.  The  distribution  of  the  damping  coefficients  over  the  body  is  also 
very  well  predicted. 


13 


Figure  8.  Comparison  of  the  damping  moment  coefficients  predicted  from  Sacks’ 
relations  with  CFD  predictions,  ANSR,  Mach  2.5,  xcg/L  =  0.5598 ,  and 
L/D  =  9. 


Figures  9-1 1  show  similar  comparisons  for  flight  velocities  of  Mach  1.8,  3.5,  and  4.5.  The 
relative  comparisons  are  very  similar  to  the  predictions  at  Mach  2.5  shown  in  figure  8.  The 
largest  Mach  number  effect  appears  to  occur  on  the  cylindrical  portion  of  the  body  rather  than  on 
the  nose.  The  relative  Mach  number  effect,  though  relatively  small,  is  captured  in  the  application 
of  the  Sacks’  relations. 


Similar  analysis  can  be  performed  for  the  force  coefficients.  Figure  12  shows  the  comparison  of 
the  individual  pitch-damping  moment  coefficients  CNq  and  CN^  predicted  using  the  form  of 

Sacks’  relations  in  equations  27  and  28  with  direct  CFD  predictions.  The  distribution  of  the 
force  coefficients  along  the  body  is  very  well  predicted  using  Sacks’  relations: 
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Figure  9.  Comparison  of  the  damping  moment  coefficients  predicted  from  Sacks’ 
relations  with  CFD  predictions,  ANSR,  Mach  1.8,  xcg/L  =  0.5598 ,  and 
L/D  =  9. 


Figure  10.  Comparison  of  the  damping  moment  coefficients  predicted  from  Sacks’ 
relations  with  CFD  predictions,  ANSR,  Mach  3.5,  xcg/L  =  0.5598 ,  and 
L/D  =  9. 
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Figure  11.  Comparison  of  the  damping  moment  coefficients  predicted  from  Sacks’ 
relations  with  CFD  predictions,  ANSR,  Mach  4.5,  xcg/L  =  0.5598 ,  and 

L/D  =  9. 


Figure  12.  Comparison  of  the  damping  force  coefficient  predicted  from  Sacks’  relations 
with  CFD  predictions,  ANSR,  Mach  2.5,  xcg/L  =  0.5598 ,  and  L/D  =  9. 
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Figure  13  shows  the  results  for  the  individual  pitch-damping  coefficients  for  various  CG 
positions  for  the  body  with  the  L/D  ratio  of  9  at  Mach  2.5  obtained  by  applying  equations  27  and 
28.  The  results  are  compared  with  CFD  results  and  with  direct  slender  body  theory  results. 


Figure  13.  Comparison  of  the  predicted  damping  moment  coefficients  with  CFD  and 
slender  body  theory  results,  ANSR,  Mach  2.5,  and  L/D  =  9. 

There  is  very  good  correlation  of  the  results  obtained  with  Sacks’  relations  and  the  CFD  results, 
and  the  results  are  a  significant  improvement  over  the  slender  body  theory. 

For  the  results  shown  in  figures  8-12,  comparisons  have  been  made  between  the  force  and 
moment  distributions  along  the  body.  This  was  done  because  these  distributions  are  easily 
extracted  from  the  computational  results,  and  they  allow  detailed  examination  of  the 
performance  of  Sacks’  relations.  It  must  be  emphasized  that  it  is  not  necessary  to  obtain  the 
pitch-damping  force  and  moment  distributions  in  order  to  apply  Sacks’  relations  because  they  are 
applicable  for  the  global  force  and  moment  coefficients  as  well.  This  is  particularly  important  to 
consider  when  the  source  of  the  pitch-damping  coefficient  data  is  from  experiment  where  only 
global  coefficients  are  typically  available.  However,  Sacks’  relations  do  require  the  normal  force 
distribution  to  determine  C2m  because  generally  only  the  normal  force  and  pitching  moment 
are  available  as  global  coefficients  from  many  sources.  In  lieu  of  more  sophisticated 
computational  methods  for  determining  the  normal  force  distribution,  fast-design  aeroprediction 
codes  (5)  should  provide  acceptable  accuracy  for  a  variety  of  flight  vehicle  geometries. 
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To  demonstrate  the  use  of  Sacks’  relations  when  only  global  forces  and  moments  are  available, 
pitch-damping  moment  coefficient  sum  data  from  range  firings  of  the  ANSR  have  been  used  to 
estimate  the  individual  pitch-damping  moment  coefficients  using  Sacks’  relations.  The  normal 
force  distribution  from  the  fast-design  aeroprediction  code  AP02  (3)  has  been  used  to  predict 
C2m  .  Figure  14  shows  the  predicted  pitch-damping  coefficients  obtained  from  Sacks’  relations 
using  experimental  measurements  of  the  pitch-damping  coefficient  sum  and  C2m  obtained 
from  AP02.  The  results  are  compared  with  direct  CFD  predictions  of  the  individual  damping 
coefficients  as  well  as  the  results  obtained  from  Sacks’  relation  using  CFD  predictions  of  the 
pitch-damping  coefficient  sum  and  C2m^  shown  previously  in  figure  10.  The  experimentally 
derived  values  of  the  individual  pitch-damping  coefficients  compare  well  with  the  predicted 
results.  The  results  also  indicate  that  the  biggest  source  of  error  is  produced  by  the  uncertainty  in 
the  experimentally  derived  pitch-damping  moment  coefficient  rather  than  the  error  in  Sacks’ 
relations.  (The  pitch-damping  data  used  here  are  from  a  highly  regarded  data  set  and  are 
representative  of  the  expected  accuracy  for  the  pitch-damping  coefficient  sum  from  range 
firings.) 


Xcg/D 


Figure  14.  Comparison  of  the  damping  moment  coefficients  derived  from  experimental 
data  or  CFD  using  Sacks’  relations  with  direct  CFD  prediction,  ANSR,  Mach 
2.5,  and  L/D  =  9. 

Both  Sacks  (1 )  and  Bryson  (2)  obtained  relations  between  the  pitch-damping  coefficients  and  the 
normal  force  coefficient  shown  previously  in  equations  1  and  2.  Relatively  speaking,  equations 
1  and  2  do  not  perform  particularly  well,  especially  when  compared  to  Sacks’  relations 
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(equations  4-9).  Figure  15  shows  a  comparison  of  the  CFD  prediction  of  the  pitch-damping 
coefficient  CNq  with  results  obtained  using  equation  2  and  Sacks’  relation  (equation  28).  As 
before,  CFD  results  were  used  to  obtain  the  aerodynamic  coefficients  required  in  equations  2  and 
28.  The  results  obtained  with  equation  2  provide  only  a  qualitative  prediction  of  the  distribution 
of  CN  compared  with  the  result  obtained  using  Sacks’  relation  (equation  28),  which  is  in  much 
better  agreement  with  the  CFD  predictions.  Instead  of  relating  the  pitch-damping  coefficients  to 
the  normal  force  coefficient  as  in  equation  2,  Sacks’  relations  relate  the  pitch-damping 
coefficient  to  each  other.  These  results  seem  to  imply  the  pitch-damping  coefficients  are  more 
closely  related  to  each  other  than  they  are  to  the  normal  force  coefficient.  Figure  16  shows  a 
comparison  of  the  pitch-damping  moment  coefficient  sum  Cm  +  Cm^  from  direct  CFD 
prediction  with  the  results  obtained  using  Bryson’s  relation  (equation  1),  which  like  equation  2, 
relates  the  damping  coefficient  to  the  normal  force  coefficient.  The  distribution  of  Cm  +  Cm^ 
is  relatively  poorly  predicted  compared  with  the  CFD  result,  although  total  moment  for  entire 
configuration  differs  by  only  25%. 


Figure  15.  Comparison  of  the  pitch-damping  force  coefficient  CN  obtained  from 

direction  CFD  predictions,  predictions  using  Sacks’  relation  (equation  28) 
and  equation  2,  ANSR,  Mach  2.5,  xcg/L  =  0.5598 ,  and  L/D  =  9. 

3.2  Flared  Projectile  Results 

The  performance  of  Sacks’  relations  was  also  examined  for  a  flared  projectile  geometry  shown  in 
figure  17.  Predictions  of  the  pitch-damping  coefficient  sum  have  been  validated  previously  for 
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Figure  16.  Comparison  of  the  pitch-damping  moment  coefficient  sum  predicted  from 
Bryson’s  relation  (equation  1)  with  CFD  predictions,  ANSR,  Mach  2.5, 
xcg/L  =  0.5598 ,  and  L/D  =  9. 


this  configuration  (14).  For  the  current  report,  additional  predictions  of  the  individual  pitch¬ 
damping  coefficients  were  also  performed  for  flight  velocities  between  Mach  2  and  Mach  5. 
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Figure  18  shows  the  predicted  deviations  of  the  moment  compared  with  the  individual  pitch¬ 
damping  coefficients  at  Mach  2.  Again,  the  deviation  over  the  body  is  small  compared  to  the 
pitch-damping  coefficients  even  on  the  afterbody  where  both  the  geometry  and  moment 
coefficients  are  changing  significantly.  Similar  results  were  found  at  the  other  flight  velocities. 


Figure  18.  Comparison  of  the  three  computed  deviations,  A, ,  A, ,  and  A3  relative  to 
the  damping  moment  coefficients,  flared  projectile,  and  Mach  2. 

Similar  to  the  results  for  the  secant  ogive/cylinder  configuration,  the  form  of  Sacks’  relations 
embodied  in  equations  25  and  26  was  applied  to  determine  the  individual  pitch-damping 
coefficients  from  the  pitch-damping  moment  coefficient  sum  and  the  pitching  second-moment 
coefficient  for  the  flared  projectile  geometry.  Comparisons  were  made  with  direct  CFD 
predictions. 

Figures  19  and  20  show  comparisons  of  the  individual  pitch-damping  moment  coefficients  at 
Mach  2  and  5,  respectively,  for  the  flare  projectile  geometry.  Very  good  agreement  between  the 
results  obtained  by  applying  Sacks’  relations  and  direct  CFD  predictions  are  found.  Additional 
results  for  the  same  flare  projectile  with  the  6°  flare  replaced  by  a  15°  flare  were  also  obtained 
and  are  shown  in  figure  2 1 .  Despite  a  nearly  doubling  of  the  pitch-damping  coefficient  due  to 
the  larger  flare,  the  Sacks’  relations  results  are  in  very  good  agreement  with  the  direct  CFD 
predictions. 
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XI D 


Figure  19.  Comparison  of  the  damping  moment  coefficient  predicted  from  Sacks’ 
relations  with  CFD  predictions,  flared  projectile,  and  Mach  2. 


Figure  20.  Comparison  of  the  damping  moment  coefficient  predicted  from  Sacks’ 
relations  with  CFD  predictions,  flared  projectile,  and  Mach  5. 
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Figure  21.  Comparison  of  the  damping  moment  coefficient  predicted  from  Sacks’ 
relations  with  CFD  predictions,  15°  flared  projectile,  and  Mach  5. 


4.  Summary 


In  conclusion,  the  results  presented  here  indicate  that  Sacks’  pitch-damping  relations  are  only 
strictly  valid  under  the  context  of  the  theory  from  which  they  were  originally  developed.  They 
do,  however,  provide  a  reasonably  good  means  of  estimating  the  pitch-damping  coefficients 
when  one  of  the  three  pitch-damping  coefficients  can  be  determined.  The  most  likely  practical 
use  of  these  relations  might  be  to  provide  estimates  of  the  individual  pitch-damping  coefficients 
using  values  of  the  pitch-damping  coefficient  sum  determined  from  some  other  source,  such  as 
experimental  data  and  engineering  estimation  approaches  or  when  the  additional  expense  of  the 
separate  CFD  computations  of  the  individual  pitch-damping  coefficients  is  not  justified.  In  some 
cases,  it  appears  that  the  error  in  applying  these  relationships  is  smaller  than  the  error  associated 
with  generating  the  initial  pitch-damping  coefficient  (such  as  with  engineering  estimation 
approaches)  from  which  the  other  two  damping  coefficients  are  derived  using  Sacks’  relations. 

Applying  Sacks’  relations  to  determine  the  individual  pitch-damping  coefficients  from  the  pitch¬ 
damping  sum  represents  only  one  possible  application  of  Sacks’  relations.  These  relations  could 
also  benefit  theoretical  developments  because  theories  for  predicting  the  pitch-damping 
coefficients  need  only  focus  on  a  single  damping  coefficient.  The  other  damping  coefficients 
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could  then  be  obtained  from  Sacks’  relations.  Such  an  approach  has  been  already  used  as  an 
estimation  procedure  for  the  damping  coefficients  ( 6 ).  In  this  work,  the  distribution  of 
CN  .  along  the  body  is  predicted  using  slender  body  theory  with  empirically  based  corrections. 
The  damping  coefficient  CNq  can  then  be  obtained  from  Sacks’  relations.  Once  the  damping 
force  distributions  are  known,  the  damping  moments  can  be  easily  obtained  by  integration  of  the 
force  loadings.  Improvements  in  the  estimates  of  CNq  from  Sacks’  relations  can  also  be 
obtained  by  correlating  the  error  A  as  well  -  an  approach  used  in  reference  (6)  to  further 
improve  the  estimates  of  CNq  . 

Finally,  in  the  currently  report,  only  axisymmetric  configurations  in  supersonic  flight  have  been 
considered.  The  theory  from  which  Sacks  derived  the  relations  considered  here  is  applicable  to 
both  wingless  and  winged  vehicles.  Further  research  is  still  required  to  assess  the  perfonnance 
of  Sacks’  relations  for  winged  vehicles  and  for  other  flight  velocity  regimes. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


Speed  of  sound 


Pitching  moment  coefficient, 


—  SrefD 


Pitching  moment  coefficient  slope  with  respect  to  angle  of  attack, 


Pitching  second  moment  coefficient  slope  with  respect  to  angle  of  attack, 


dC2m 


Pitch-damping  moment  coefficient  slope, 


Pitch-damping  moment  coefficient  slope, 


Cmq  +  C  Pitch-damping  moment  coefficient 

C„  Side  moment  coefficient 


Normal  force  coefficient, 


Sref 


Normal  force  coefficient  slope  with  respect  to  angle  of  attack,  - — 

da 


Pitch-damping  force  coefficient  slope, 


Pitch-damping  force  slope, 


CNq  +  CN(.  Pitch-damping  force  coefficient 
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D 

e 

e,f,g 

F 

H 

Hc 

L 

M 

M 

P 

q 

qT 

r 

Re 

Ro 


S 


ref 


s 

Sc 

u,v,  w 
U,V,W 


V 


Xe ,  Ye ,  Ze 


x 

x 


X 


eg 


Reference  diameter 

Total  energy  per  unit  volume 

Flux  vectors  in  transformed  coordinates 

Force 

Source  term  in  Navier-Stokes  equations  due  to  rotating  coordinate  frame 
Source  term  in  Navier-Stokes  equations  due  to  cylindrical  coordinate  formulation 

Body  length 
Mach  number 
Moment 
Pressure 

Transverse  angular  rate  of  body  as  used  in  flight  mechanics  equations 
Transposed  vector  of  dependent  variables  as  used  in  Navier-Stokes  equations 

Radial  coordinate 
Reynolds  number,  a^/^D///, 

Helix  radius 


Reference  area, 


S 


ref 


STD2 

4 


Viscous  flux  vector 

Viscous  terms  due  to  cylindrical  coordinate  formulation 
Velocity  components  in  the  x,^,r  directions 
Contravariant  velocity  components 

Freestream  velocity 
Earth-fixed  coordinates 

Axial  location  along  body  from  nose 

Integration  variable  associated  with  axial  location  along  body  from  nose 
Axial  location  of  center  of  gravity  from  nose 
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Greek  Symbols 


a  Angle  of  attack 

a  Angular  rate  associated  with  angle  of  attack 

P  Yaw  angle 

y  Cosine  of  total  angle  of  attack 

8  Sine  of  total  angle  of  attack 

A  Deviation  or  error  in  Sacks’  relations  for  force 

A  Deviation  or  error  in  Sacks’  relations  for  moment 

//  Viscosity 

E,,  r),  8,  Transformed  coordinates  in  the  Navier-Stokes  equations 

p  Density 

(f>  Circumferential  coordinate 

<x>  Angular  velocity  about  longitudinal  axis 

Q  Angular  rate  associated  with  coning  and  helical  motions 

Subscripts 

go  Quantity  evaluated  at  freestream  conditions 
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